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We present a general relativistic model of a spherical shell of matter with a perfect fluid on 
its surface coupled to an internal oscillator, which generalizes a model recently introduced by the 
authors to construct a self-gravitating interferometer [J| . The internal oscillator evolution is defined 
with respect to the local proper time of the shell, allowing the oscillator to serve as a local clock 
that ticks differently depending on the shell’s position and momentum. A Hamiltonian reduction is 
performed on the system, and an approximate quantum description is given to the reduced phase 
space. If we focus only on the external dynamics, we must trace out the clock degree of freedom, and 
this results in a form of intrinsic decoherence that shares some features with a proposed “universal” 
decoherence mechanism attributed to gravitational time dilation Q]. We note that the proposed 
decoherence remains present in the (gravity-free) limit of fiat spacetime, emphasizing that the effect 
can be attributed entirely to proMr time differences, and thus is not necessarily related to gravity. 
Whereas the effect described in [2| vanishes in the absence of an external gravitational field, our 
approach bootstraps the gravitational contribution to the time dilation decoherence by including 
self-interaction, yielding a fundamentally gravitational intrinsic decoherence effect. 


I. INTRODUCTION 

Time dilation is one of the most profound consequences 
of relativity theory. In classical systems, time dilation ef¬ 
fects are fairly well understood, but in quantum systems 
there are still many questions that remain unanswered. 
One such question is how to properly incorporate the ef¬ 
fects of time dilation into the quantum evolution of com¬ 
posite systems, either in the limit of flat spacetime or in 
situations where gravitational effects are significant. 

In a series of papers by Pikovski et al. i-i , for in¬ 
stance, a “universal” decoherence mechanism was pro¬ 
posed for composite general relativistic systems, due to 
gravitational time dilation. In Q, Pikovski et al. present 
an approximate quantum description of such a compos¬ 
ite system that as a whole behaves as a point particle 
(located at the center-of-mass of the system) with a well- 
defined proper time. The system is placed in the grav¬ 
itational field of the earth, and has internal degrees of 
freedom that are defined in the rest frame of the system. 
Expressing the proper time derivative in terms of the lab- 
frame time t induces a coupling between the internal de¬ 
grees of freedom and the center-of-mass coordinate, and 
if one only keeps track of the center-of-mass dynamics, 
tracing out over the internal degrees of freedom leads to a 
novel form of the “third-party decoherence” described by 
Stamp Q. It is this effect, and variations on the theme, 
that we focus on in this paper. 

Rather than make use of the same model used by 
Pikovski et ah, we will explore the same ideas with a 
model that generalizes the self-gravitating spherical per¬ 
fect fluid shell introduced recently by the authors [l|. 
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Whereas our original model was introduced to study the 
consequences of general relativity on massive interferom¬ 
eters, here we extend the model to include an “internal” 
harmonic oscillator, to analyze the quantum structure of 
composite relativistic systems. By “internal,” we mean 
that the harmonic oscillator is described by an internal 
coordinate q that oscillates in an abstract space that is 
not part of the spacetime; such an internal degree of free¬ 
dom could represent the values of a single spherically 
symmetric mode of an oscillating field confined to the 
surface of the shell, for instance. The internal coordinate 
oscillates harmonically with respect to the proper time of 
the (external) shell position, and therefore the oscillator 
serves as a clock, evolving based on the local flow of time 
determined by the external motion. 

We should keep in mind that the internal oscillator 
contributes to the external shell dynamics as well, which 
in turn affects the spacetime; in other words, the very 
ticking of our clock influences the manner in which it 
ticks. This is especially relevant in the quantized sys¬ 
tem, because uncertainties in clock readings become inti¬ 
mately connected with uncertainties in spacetime geom¬ 
etry. 

We will explore some of the ambiguities associated with 
the quantum theory of this generalized shell system in re¬ 
duced phase space, and then relate an approximate form 
of our reduced Hamiltonian with the Hamiltonian pre¬ 
sented in Q . We exploit this parallel to demonstrate time 
dilation decoherence in our system. Even in the (gravity- 
free) limit of flat spacetime, we observe that when the 
fluid pressure is nonzero, decoherence (in the position ba¬ 
sis) remains present, because of the acceleration caused 
by the pressure. This is a reflection of the fact that the ef¬ 
fect proposed by Pikovski et al. results from proper time 
differences alone, and as such is not necessarily related 
to gravity. 

In the general case, our shell model includes gravita- 
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tional self-interaction corrections to the time dilation de¬ 
coherence, such that the decoherence is altered by the 
manner in which our shell and its clock influence the 
state of their own geometry. We find that even with¬ 
out pressure, the self-gravitation of the shell leads to the 
nonzero acceleration required to produce the position- 
basis time dilation decoherence. We interpret this “self¬ 
decoherence” as a fundamentally gravitational effect. 


II. CLASSICAL ACTION 


For context, before adding an internal oscillator, the 
self-gravitating spherical perfect fluid shell model intro¬ 
duced in [I| is described by the action 1^;+ Iq, where Ix 
is the shell action 


with all quantities evaluated on the shell history, and Iq 
is the Einstein-Hilbert action 

1 = ^ j (2) 

Here superscripts on the metric determinant g and the 
Ricci scalar TZ indicate that these are constructed from 
the full (3-|-l)-dimensional spacetime metric components 
We make use of the ADM form of the metric in 
spherical symmetry, 

g^^dx^dx'^ = -N^df + {dr + N^dtf + R^dn"^, (3) 


where N is the lapse function, N'^ is the radial com¬ 
ponent of the shift vector, and and are the only 
nontrivial components of the spatial metric @ . It is then 
clear that R is the “radius” of the shell, obtained from 
the area 47ri?^ of symmetry two-spheres. The shell con¬ 
tribution Ix is analogous to a free relativistic particle ac¬ 
tion, except with a mass M that depends on the position- 
dependent metric function R; the function M{R) serves 
to parametrize the relationship between the density cr = 
M{R)/I'kR^ and pressure = —M'{R)/8nR of the 
fluid. 

We add an internal oscillator to our shell with the ac¬ 
tion 


I,= 




( 4 ) 


with T being the proper time evaluated on the shell his¬ 
tory, and q being an internal coordinate that does not 
take values in the (external) spacetime. The quantity k 
is related to wq, the natural frequency of the oscillator, 
via k = mujQ. The action Q is manifestly invariant un¬ 
der coordinate transformations, as it only makes use of 
the proper time of the shell. 

It simplifies the description to parametrize the shell 
history with the coordinate time t, such that the classical 


shell motion is defined by a trajectory r = X{t). We can 
use the shell 4-velocity = {dt/dT){l, A, 0,0) to express 
the proper time differentials as 

dr = —u^dx^ = rdt (5) 

and 

- = = r (6) 

Now q is being treated as a function solely of coordinate 
time t, to reflect our choice of parametrization. The 4- 
velocity normalization = — 1 then implies that one 
can express the derivative of the proper time with respect 
to the coordinate time as 


T = dr dm - m{m P XY5{r - X) 


= \jm-m{m + xY. 


( 7 ) 


An overhat denotes that a quantity is to be evaluated on 
the shell history; likewise, it is understood that the over¬ 
dots denote coordinate-time derivatives along the shell 
trajectory. 

If we define the original shell Lagrangian as 


/' 


Cx = - I dr d m - L'^{m + XfM{R)5{r - X) 


= -fM 


( 8 ) 


and the oscillator Lagrangian as 


1 
1 

2 


dr T 


m ( ^ ) — kq"^ 


5{r - X) 


m— - kTq‘ 

f 


then the shell-oscillator action is given by 

I shell = JdtC = J dt {Cx+Cq). 


(9) 


( 10 ) 


III. HAMILTONIANIZATION 


We can Hamiltonianize the shell-oscillator system with 
the Legendre transformation R — PX + pq — C. The 
momentum p conjugate to the internal coordinate q is 
given by 

p = ^ = m [ dr ^S{r — X) = m^, (11) 

dq J T T 


and that the momentum conjugate to the shell position 
X is given by 


P = 


dC - dr 
—^ = -M^ 
dX dX 


-m 


M ■ 


dr 

df / - 

-^ [M 

dX V 


1,2 

-h -kq^ 

2m 2 ^ 


H, 


-)■ 


dr 


hq^ — 

2^® dX 


( 12 ) 
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where the internal “clock” Hamiltonian Hq = /2m + 

kq^/2 takes the form of a (free) harmonic oscillator. In¬ 
troducing the notation M = M + Hq, one can observe 
that our shell-oscillator system becomes very similar to 
the shell system without the oscillator, subject to the 
transformation M —>■ M. 

More explicitly, the shell momentum P can be ex¬ 
pressed as 


functions (L and R). Accordingly, one can obtain an un¬ 
constrained description of the system by making a coordi¬ 
nate choice, solving the constraints for the corresponding 
gravitational momenta, and inserting the solutions into 
the Liouville form P on the full phase space, 

P = PSX+pSq+ ( dr {ttl 6L + TTfiSR). (19) 


P = dr 




M 


-h a) 


--5{r-X), (13) 


from which we can solve for X to obtain 
NP 


X = dr 


NP 


Wp'^ + l^m^ 


-N’-j Sir-X) 

- iV”. (14) 


The Legendre transformation then gives us the Hamilto¬ 
nian 


n = PX +pq- C = J dr {NH/ + N^H/), (15) 

with the definitions 

H/ = VL-2P^+M^6{r-X), 

H/ = -PSir-X). (16) 

We remind the reader that M = M + Hq, so our clock 
Hamiltonian adds to the (position-dependent) shell mass 
M to alter the Hamiltonian constraint from the form it 
took in [I|. 

Hamiltonianizing the gravitational sector as well leads 
to the total action 

I = J dt (^PX + pdj 

-I- J dtdr nR + ttlL — NHt — N^Hr'^,(17) 

for Ht = H/ + H^ and Hr = H/ + H//, such that 

rjG ^ AttI TTLTTfl / RR'V {R'f L 

* 2i?2 R \ L j 2L 2’ 

H° = R'TTR-LTrl (18) 

IV. REDUCED PHASE SPACE QUANTIZATION 

Since we are working in spherical symmetry, the metric 
itself has no actual degrees of freedom, because although 
there are only two gravitational constraints {Ht = 0 and 
Hr = 0), there are also only two independent metric 


This amounts to a pullback of the full Liouville form 
to the representative hypersurface defined by the coor¬ 
dinate choice. From the Liouville pullback, denoted by 
P, we can deduce the canonical structure of the reduced 
phase space, which only depends on the shell-oscillator 
variables X and q (and their momenta). 

To solve the gravitational constraints, first consider 
the following linear combination of the constraints, away 
from the shell: 


for 


- —Ht - —Hr =M', 

L RL 


M{r) 


A , R R{R'? 

2R 2 2L2 


( 20 ) 


( 21 ) 


The quantity A4{r) corresponds to the ADM mass H 
when evaluated outside of the shell, and vanishes inside 
the shell. We can now solve for the gravitational mo¬ 
menta TTi, TTft away from the shell. The result is 


Assuming a continuous metric and singularity-free 
gravitational momenta, we can integrate the gravita¬ 
tional constraints {Ht = 0 and Hr = 0) across the shell, 
from which we obtain the jump conditions 

, V P 

AR' = Attl = -y, (23) 

Jri L 

where V = VP'^ + M^. We use A to denote the jump of 
a quantity across the shell (at r = X{t)). 

The coordinates we will use resemble the Painleve- 
Gullstrand coordinates {L = 1,R — r}, though the 
jump conditions force us to include a deformation re¬ 
gion {X — e < r < X) near the shell. By inspection, the 
required metric function R can be generalized as 

R{r,t)=r-j^Vg(^^^y (24) 

for a function g having the properties 


lim = 1 

2 -S- 0 + dz 

(25) 

lim =0, 

2 ^ 0 - dz 

(26) 



















4 


from which it follows that 



lim R' (X — e) 
£->■0 

-4 

(27) 

lim R'{X + e) 

= 1. 

(28) 


By inserting the gravitational momentum solutions 
(1^^ associated with the coodinate choice (|24l) into the 
jump equations (|23l) and squaring, one finds 

_ iQ- 2 i OTJ 

H — (29) 

This implies that the unreduced momentum P is im¬ 
plicitly defined as a function of the reduced phase space 
quantities X, q, p, and H: one can easily solve (l2^ to 
obtain 



The ± in (150)1 indicates whether the shell is outgoing (-(-) 
or ingoing (-), with respect to our choice of coordinates. 

Let us now calculate the pullback of the full Liouville 
form to the representative hypersurface defined by our 
coordinate choice. The condition L = 1 and the fact that 
R = r outside of the deformation region implies that 

P = P5X+pdq+ f driTjiSR. (31) 

Jx-e 

We can then simplify the remaining integral by changing 
the integration variable from r to v = R', which yields 

r^'- (i-v) 

/ dr tthSR = XdX / dv = + 0{e), (32) 
Jx-e Jl - 1 

with R'_ being R' evaluated just inside the shell. Now 
the integration is trivial, and we can easily obtain the 
desired Liouville form pullback, 

P = PcSX + p5q, (33) 


It can also be shown [I| that the expression (l34t for 
the reduced canonical shell momentum is equivalent 
to the form analogous to the expression given by Kraus 
and Wilczek Q, 


Pc = -VWX - X\a. 


' X + V - P -V2HX\ 


(35) 


despite the different minus sign placement. 

To gain some intuition for how the presence of the oscil¬ 
lator alters the reduced dynamics, let us consider the flat 
spacetime limit (X —>■ oo) of the system defined by (1351) . 
Keeping in mind the similarities with the relativistic- 
particle-like structure of our shell system, it should be 
unsurprising that in this limit (1351) becomes 


Pc = ±\/p2_M2 = ±y'p2_ (36) 

and therefore the Hamiltonian is given by 

H = ^P2 + M2 = (37) 


In the nonrelativistic regime (i.e. small Pc), the Hamil¬ 
tonian can then be expressed as 


H 


M H„ ~\- 


{m- 

,2 




= M + l-kq^ + ^ 


2m 


+ 


PJ 


2 (m + 


.(38) 


The last term in this approximate Hamiltonian is an ef¬ 
fective coupling between the internal oscillator variables 
{q and p) and the external shell variables {X and Pc). 
The coupling is of course produced by the fact that the 
internal “clock” oscillates harmonically with the shell’s 
proper time, the flow of which is influenced by the exter¬ 
nal variables. 

Even in the flat spacetime nonrelativistic limit, one can 
tell from the appearance of the clock Hamiltonian Hq in 
the denominator of the last term in (1381) that exact quan¬ 
tization will require nonstandard techniques. The Hamil¬ 
tonian (1551) leads to the following Schrddinger equation, 
in the coordinate basis: 


with the reduced canonical momentum for the shell po¬ 
sition satisfying 



M + -kq^ 


— 

2m dq'^ J 


Pc = - V2HX + X\n 



V + P 
X 



(34) 


1 d 
2dX 


1 

M + \kq^ 


1 

2m dq^ _ 


dX 


4-. (39) 


This expression gives an implicit definition of the Hamil¬ 
tonian H on the reduced phase space, as a function of the 
shell-oscillator variables {X and g), along with the mo¬ 
menta that are conjugate to them in the reduced phase 
space (Pc and p, respectively). 


The factor-ordering in the last term of (1391) was chosen to 
make the differential operator Hermitian, but there is still 
some ambiguity in the meaning of the bracketed factor 
between the X-derivatives, since the formal expression 
has g-derivatives in the denominator. 
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Instead of working directly with the Hamiltonian (l38t , 
we will work with an approximate Hamiltonian that al¬ 
lows for a simpler analysis. Before doing so, however, let 
us consider the lowest-order corrections to (1551) due to 
gravitational self-interaction. This can be accomplished 
by inserting the ansatz H = Fo{X)+Fi{X)Pc+F 2 {X)P'^ 
into (|35)l and expanding for large X and small Pc- If 
we assume that the X-dependence of M does not sig¬ 
nificantly affect the large-X behaviour of the functions 
{Fj(X)}, then we find Q 

A>2 

= M-^+0{l/X^) 

F, = + (40) 


exploit the similarity enough to demonstrate an analo¬ 
gous decoherence in our system, as will become clear in 
what follows. 

Represented by a density operator, the full state p 
obeys the von Neumann equation, 

p = -i [H, p]. (44) 


We can then change the frame to primed coordi¬ 
nates, as in 0, which are defined by p'{t) = 

h(X,P,) = Trg[H,gPg{0)]. 
We are assuming that the initial state of the system is of 
the product form p(0) = px{0)pq{0), i.e. initially uncor¬ 
related. Denoting the average clock energy Tr, [F[qPq{0)] 
by Eq and the shell part of the interaction by 


r 


p2 

c 

2M2 


-b 


M 1 p 


(45) 


We will now assume that the shell mass M is suffi¬ 
ciently larger than the clock energy F[q that we can ex¬ 
press M as M ^I -b F[q/M^ and expand the Hamiltonian 

in powers of Hq/M. To lowest order, the system then de¬ 
composes into a more standard form, given by 

H = Ho + H^q=H^+Hq + H^q, (41) 

with the approximate shell Hamiltonian 


we obtain the expression h = T(X, Pc)Eq. In equation 
(SSI), as well as equations SI and (H51) . the bracketed 
term originates from self-gravitation. The transformed 
von Neumann equation is 

p\t) = imt)+h'it),p'{t)]-t[H',it)+H',^it),p'it)] 

= (46) 

with h'{t) = h{X'{t),P^{t)). If we integrate and iter¬ 
ate equation (l46l) . we are led to the integro-differential 
equation 


H^=M 


Pc 

1 

2 

h-^Pc- 


2M 

3X 

3\ 

X 

18X 


(42) 


the internal “clock” Hamiltonian Flq = p^/2to -b kq^/2, 
and the interaction 


p2 

TT _ C TT 


^xq 


2M2 


-—Hq - 

9X « 3mV ^ 


(43) 


which is induced by the clock oscillation being defined 
with respect to the proper time of the shell. 


V. TIME DILATION DECOHERENCE 

The decomposition SB) is of the same form as the 
one used recently by Pikovski et al. Q to demonstrate 
decoherence due to gravitational time dilation for com¬ 
posite systems, though the interpretation of the system 
variables is different. The main difference is that the 
“external” coordinate of our system is the shell radius 
instead of the (somewhat ill-defined) center-of-mass co¬ 
ordinate. In our case, however, the self-gravitation of 
the shell-plus-clock system is taken into account, so both 
the shell and the clock influence the spacetime geometry, 
which is therefore no longer fixed. Nonetheless, we can 


P\t) 


-i [K^{t)-h'{t),p'{0)] 




Hs),p'{s) , 


(47) 


using the definition h{t) = — h'{t). At this point 

Pikovski et al. trace over the internal variables, which for 
us describe the clock, and make use of the Born part of 
the Born-Markov approximation, keeping only terms up 
to second order in the interaction Hamiltonian H^q, and 
replacing the p'{s) in the integral by p(,(s)p(^(0). This ap¬ 
plication of the Born approximation assumes weak cou¬ 
pling, but in contrast to the full Born-Markov it does 
not ignore memory effects. For a detailed discussion of 
this approximation, see @ . The reduced equation for the 
shell variables is then given by 


P'x(t) =Trq [p'{t)] 

Hs),p'{s) I 

= - ds Trq {{Hq - Eqf [P'(t), [r'(s), p'(s)]]} 

= -{XEqf f ds[T'{t),[T'{s),p',{s)]], (48) 

Jo 

with the notation r'(s) = r(-^^(s)) ^’c('S)) and 


i^Eqf 


Tr„ 


{Hq 


EqfPq ( 0 ) 


(49) 
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One can then transform back to the unprimed frame, whereby the substitution s—>■ t—s leads to the expression 


J 


Pxit) = -i [Ha; +TEg,pa;{t)] - {AEqf 


In general, the reduced evolution equation (EOI) ex¬ 
hibits decoherence due to the nonunitary contribution of 
the last term on the right. Under some special circum¬ 
stances this term vanishes, leaving the reduced system 
to evolve unitarily; for example, such a circumstance oc¬ 
curs for initial states that are eigenstates of the internal 
(clock) Hamiltonian, of course then the oscillator is not 
much of a clock, as it never changes with time (modulo 
a phase). 

VI. DISCUSSION 

In the last section, we demonstrated intrinsic decoher¬ 
ence due to time dilation. The decoherence basis is in 
general a combination of position and momentum, re¬ 
ducing to purely position in the limit of negligible shell 
momentum. In the system described by Pikovski et al. 
0, which includes spacetime curvature caused by the 
external gravitational field of the earth, the part of the 
time dilation decoherence that involves the position basis 
should vanish in the absence of the earth’s gravitational 
influence: without the earth, the center-of-mass coor¬ 
dinate they use defines the origin of an inertial frame, 
and in that frame the proper time associated with the 
center-of-mass coordinate is equal to the coordinate time 
However, in our system, this decoherence is present 
even in the limit of flat spacetime (ignoring both exter¬ 
nal gravitational fields and self-gravitation), because of 
the nonzero acceleration of the shell due to the position 
dependence of the mass (M = M{X)). We then see 
a confirmation that the time dilation decoherence pro¬ 
posed in 0 is not necessarily related to gravity, but pro¬ 
duced by proper time differences in composite systems 
with nonzero accelerations. 

Although a simpler form of the shell model we employ 
was previously used by the authors to examine the pos¬ 
sible origins of the typ e of gravitational decoherence ex¬ 
plored by Penrose [ij and Diosi m , the self-gravitation 
described here results from classical general relativity 
alone, and is thus not related to speculative models that 
produce self-gravitation by altering quantum mechanics. 
It is still an open question whether or not Penrose-type 
gravitational decoherence can be demonstrated within 


f ds [P, [P, pait - s)] . (50) 

0 


canonical quantum gravity, without introducing any new 
physics (for further discussion, see i)- 

Another observation we can make is that even in 
the constant M limit, where the fluid pressure van¬ 
ishes, equation (l50l) indicates that the time dilation de¬ 
coherence remains present, in this case because the self¬ 
gravitation produces a nonzero acceleration of the shell 
position. Since the decohering term in the reduced evo¬ 
lution equation (1501) varies with the square of the clock’s 
energy uncertainty, increasing the energy uncertainty of 
the clock enhances the decoherence. In our system the 
clock energy contributes to the ADM energy of the space- 
time, so uncertainties in the clock energy contribute to 
uncertainties in the spacetime geometry, which in turn 
lead to the type of self-decoherence mentioned in the in¬ 
troduction. Conceptually, such an effect should occur 
for any composite general relativistic system that has in¬ 
ternal motion that (classically) evolves according to the 
proper time associated with the system’s external mo¬ 
tion, since the alteration of the local flow of time caused 
by the system’s influence on its own spacetime geometry 
induces an effective coupling between the internal and 
external degrees of freedom. It is this effect that is fun¬ 
damentally gravitational in nature, as it is present even 
in the absence of any other interactions. 

We have therefore arrived at a type of intrinsic de¬ 
coherence similar to the “third-party” decoherence de¬ 
scribed by Stamp Q, though in contrast to the use of 
the earth as the third party as proposed by Pikovski et 
al. 0, we have bootstrapped the idea by incorporating 
gravitational self-interaction, effectively producing third- 
party decoherence without the third party. 
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